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The complete phase diagram of Random Energy Model (REM) is obtained for complex tempera- 
tures using the method proposed by Derrida. We find the density of zeroes for statistical sum. Then 
the method is applied to Generalized Random Energy Model (GREM). This allows us to propose 
new analytical method for investigating zeroes of statistical sum for finite-dimensional systems. 
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Random Energy Model (REM) [1] is one of the most popular contemporary models of statistical mechanics. Besides 
H ■ its direct applications in spin-glasses the model has a wide range of employments in very different areas of the modern 
theoretical physics and biophysics. Recently, several new, and more intimate applications have been founded in the 
theories of well-developed turbulence [2], and strings [3-5]. On the other hand, it is well-known that the model allows 
to consider the fundamental problems of information theory in the language of statistical mechanics. In particular, 
problems connected with transmission of information through a noisy channel can be formulated and effectively solved 
in the language of REM [6,7]. This astonishing range of different applications cannot be accidental one: It seems, 
that REM is a paradigm for a complex system physics, or at least contains all essential ingredients of that physics. 
In other words, the nature sings his sacral songs using REM's language. 

The thermodynamical structure of REM has simple structure but at the same contains all essential ingradients 
of phase transitions [El Ojl^ ]: There is well defined critical temperature T c which can be obtained from clear and 
semi- intuitive physical reasons. For T < T c the system is frozen in the spin-glass phase with one level of replica 
symmetry breaking. Besides this rough thermodynamical phase transition there are some " mild" , mesoscopic phase 
q transitions in the domain (T,T C ), T > T c p|Jl2[. The temperature T strongly depends from the concrete form of 
energy's distribution; it can be infinite for the gaussian distribution [jjj but it has some finite value for REM generated 
by more physical dilute hamiltonians p2[ . This rich and interesting spectra of phase transitions waits its complete 
and exhaustive investigation, and we hope that the present work will shed some light to this question also. 

It is well-known that structure and properties of phase transitions can be investigated through analytical properties 
of thermodynamical quantities in the complex plane of temperature or/and magnetic field [Q Q Jl6| |. Indeed, if 
a phase transition is assosiated with singular behavior of thermodynamic potential in the thermodynamic limit (in 
• our case it is free energy or statistical sum), then by consideration of its analytical properties in the complex plane 
the important physical information can be obtained. The method was proposed by Yang and Lee |l4j (for the case 
of complex magnetic fields), and Fisher Jy| (for the case of complex temperatures). It has large variety of different 
' applications in statistical physics; furthermore it is one one of the most exact and powerful methods for estimating 
critical indices and constructing phase diagrams. In particular, the method is applied for investigation of phase 
transitions in disordered or strongly frustrated systems because sometimes all other methods are not very informative 
or hard for realization. Resently with help of this method a dense domain of phase transitions has been obtained in 
I ■ a non-disordered but fully-frustrated model [|l6| . 

REM has been investigated in complex field/temperature plane [8,10]. The goal of the present paper is to continue 
this program, and to derive the complete phase diagram. Furthermore, we propose a new method for investigating 
analytical properties of thermodynamical quantities of finite-dimensional systems. There is another, perhaps more 
fundamental motivation for the present investigation: it was obtained resently that REM at complex temperatures is 
closely connected with strings M Q . 

In disordered systems there are averaging with the canonical Gibbs distribution at given temperature, and averaging 
by frozen disorder. In REM and related models energy levels are considered as random quantities, and the distribution 
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is derived from a concrete " microscopic" Hamiltonian jl| || . 
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Here Si are ±1 spins, are quenched couplings with normal distribution and a proper scale 



P(3h..i p ) = Sl^exp^CZjL.iJN, (2) 
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where = N\/p\(N — p)\. It is possible to prove, that in the limit of large p all energy levels are distributed 
independently with probability 



Let us define partition function and free energy. Consideration of complex temperatures is some tool for investigating 
analitical properties of thermodynamical quantities. In particular, the usual definitions of statistical physics should 
be adapted to this purpose. Partition function is defined as usually: 
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Further, at the thermodynamical limit different expressions like 

< In | Z |>,< In | ReZ |> (4) 

are identical. In a slightly cavalier fashion we will refer to the quantity (In |ReZ|) as free energy, taking into account 
that for the case of real temperatures it is related simply to the usual free energy. 
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(ln|ReZ|)= / dt- (/ dt- )= / dtlnte- t + I lnid( e -*l RcZ l) (5) 



where A = flvN, and (...) means averaging by the normal distribution. Further, 

(e-*l ReZ l) - (e-*) = 4= / dze- z2 / 4 f (zVt, A) 2 
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and 
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here we have ft = Pi + fo, A = Ai + i\ 2 = VN (3 

We see that \f(A, A)| is always less than 1. When \nA^> the f(A, A) is exponentially suppressed and we get 
f(A, A) « 1, while In A <C U. To find asymptotics it is convenient to use another representation, like [8] 



f{A X) = ±l dxe^r- I dye^y-^ Xv+ ^= (8) 
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In the formula (4) the integral is taken along the line, which passes the point x = on the right side. To calculate 
asymptotic we should move the integration line at left to catch the saddle point. There is a simple pole at point 0. It 
is easy to check, that the corresponding residue equals 1. 

/ico — s 
dxe—* — xlnA g{x,A) (9) 
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where 
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g(x,A)= JyyX-leiycosXi/X^ny/A ( 1Q ) 

Jo 

This representation is useful, when the saddle point x belongs to the interval (0,-1). Function g(x, A) has a compli- 
cated singularity at point -1. We need in asymptotic at A <~ VA^hi A <~ N: 

/(A,A)-l«exp[-5^] (11) 
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where we missed pre-exponents. 

There are another 2 modes for f(A, A), besides (7). In order to calculate them we arc considering derivatives: 
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Eventually we have (after shift of integration line) 
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The similar approach gives 



/ A « 1 + z^(e A2 / 4 + e A2 / 4 ) « 1 + zAe^-^l* (14) 



fk^-^\ dxe^-^ A+ ^g(x,A) (15) 
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So we have three phases 
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rfz exp (-^ 2 /4 + 2 Jv ^e- L r- a ) = exp (-{ 2 4 w e-r- ) (16) 
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(17) 
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We obtain the following expressions for free energy. For the paramagnetic (PM) phase, without any magnetization, 
zeros of partition sum: 

(ln|RcZ|) =iVln2 + ^l^tl (19) 
For the pin glass (SG) phase with some replica symmetry breaking: 

(ln|ReZ|) = NVhP2l3i (20) 
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For the phase without any magnetization, but with some zeros for partition sum (further abbreviated as LYF; Lee- Yang 
@, and Fisher fTf]): 

(\n\RcZ\) = -Nhi2 + ^-N (21) 

Let us consider now the diluted version of the model. If we consider again our hamiltonian (|l|), only with ±1 couplings, 
then chose randomly among the all possible sets of indices i\..i v some aN 

aN 

H = - jix-ip^-'i, ( 22 ) 

(ii<ia..<ip)=l 

Now instead of the normal distributions we have 
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p(E) = — I dEie -EE 1 -aN\n.ahEi (33) 
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Here a > 1 is a ratio of couplings number to spins number. For function f(A, A) we have 



2^" J — z'oc */ — OG 

-i />ZOO />oo 1 , i 0l 1 _ i @2 

— I dxe aN In cosh f3x+x\n A / e |(j, ^ +y ~ W )y- x dy = 

2tt J- loo Jo 



1 

2^ 



Again pole at point gives 



This asymptotics is correct, when 



To calculate another types of asymptotics we have to consider the derivatives f' A , f'j^. 
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(24) 



f(A, A) « 1 - e^^^^-^^tanh/ta = -^4, (25) 
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tanh0>^— (26) 
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We have to move integration line left to catch the sadle point. Again we intersect the pole at point 0. 

i A 

f (A A) Rj 1 + — — (e aN ln cosl1 @ + e aN ln cosh @ ) 
This asymptotics is correct if for a solution 
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we have the constraint Rex < 0. 

Let us consider now the third type of asymptotics. The consideration of second derivative reads: 

fj_ = L dEe -EE 1 p+ a N\nco S hf3E 1 + !f(e> 3E +e' 3E )( e 2l3E + + 2e 2 ^ E ) = 

= -— / e aNlncosh ^ +2 ^- E ^ nA g(x)dx 

-i pioo 

- — / e aNlncosh ^ +2 ^- EllnA g(x)dx 
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The shift of integration line intersects a pole, and gives 
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where 
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Again we have 3 mods. 
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c -0 = 1 / e _ z 2/ 4+J!i 2JVt[ e «Wlncoah(/9) +e aiVlncoB(/3) ] 
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e -^ = -!- / exp [-z 2 /4 - 2 w e ] 
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e~+ = —7= / cxp (-z 2 /4 - z 2 2 4 A^ 2 e aAr ln cosh 2ft ) 

Eventually we have three expressions for free energy. For the para magnetic phase: 

< ln|i?ez| >= Aln2 + aiVi?elncosh/3 

For the LYF phase: 

N aN 

< \n\Rez\ >= — In 2 + — lncosh2/3i 



For the SG phase: 

< \nRe\z\ >= aNfii tanh/3 c (39) 

where for the j3 c we have 

a[P c tanh p c - In cosh C ] = In 2 (40) 

The search of other possible phases like 

aN N 

In cosh 401 H In 2 (41) 
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give s the same equations, as for the case of real T. So we have not any new phase, besides those obtained in Eqs. 
(p7H39| ). The boundary lines between 2 phases one can be found from the coincidence of 2 equations . 
The boundary between LYF and SG phase is the line 

/3 1 =/3 ,/3 </5 2 <oo (42) 

where the multicriticity point (/?i,/?o) lives on the intersection of line (E3) with the boundary SG, LYF. 
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sin p 2 = n (43) 

1 - tanh 2 Pi 

For the case of ferromagnetic couplings we have the ferromagnetic phase with the old expression with Pi ^ 0, fa = 0. 
Let us consider the density of partition zeros near the Pi axes. For the density of zeros we have [9]expresion 
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For the LY phase we have 
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Near the PM-SG transition line we obtain: 



7^(^52+ 7^2) <^\Z\> (44) 



p(Pl,02) = - —CnTvi ( 45 ) 



p{Pi,P2) = - , a fR , 2 S(p 1 +p 2 -p c )(p 1 -p c ) (46) 
tt cosh(pi) z 

Now we consider the phase structure of Generalized Random Energy Model in the domain of complex T . The SG 
phase transition conditions are the old ones, with the Pi instead of l/T. 

The conditions for phase transition to 3-rd (LYF) phase is more complicated. We have to consider an expression 
like 

< nnexp^[ e A ^+ A <^ +6*-*-+^] > (47) 

while expanding the exponent (to get LYF phase), we found 

< II II 1 A 2 /4[e^ +J ^+^ +J ^ E ^} > (48) 

a 0a 

This expression showes, that if at some level of hierarchy exist LY phase, then higher could not exist PM, so stays 
only the SG one. 

Let us consider the diluted version of GREM [13], with infinite number of hierarchy M . We can consider the case 
of large M with smooth distribution of zj. (number of couplings at level k) and (number of k-th level branches is 
2^*=). In this case we can introduce continuous variable v = jj between and 1, labeling the level of hierarchy and 
define distributions 

z k = dz = zdv, N k = dN = n'(v)dv dv = — (49) 
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where n(v) is a given function (entropy in bits-s). The variable v ( < v < 1) parameterizes the level of the hierarchical 
tree and z is just a parameter (for our spin system z is a total number of couplings and parameter v labeling the level 
of hierarchy) . 

Similarly to the case of dilute GREM at real T we found: 
0F 

- t jj- = z(l - v 2 (0))Re\ncosh0 + n{v 2 (0))ln2 + n(v 2 (0))ln2/2 (50) 

vi 08) 



+z{v 2 (0)-v 1 (3))/2\ncosh20 1 +z0 1 / dv g{^--) 
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where v 2 {0) 1 v\{0) are defined from the equations 

2zi?elncosh/3 + n> 2 (/3))ln2 = zlncosh2/3i (51) 

zi?elncosh/3 + n'(wi(/3))ln2 = 1 zg(z/n'(v 1 )) 

and function g(x)hom the 

1(1 + g) ln(l + g) + 1(1 - g) ln(l - g)\ = ^ (52) 

For the case of Edwards- Anderson model placed on d-dimensional hypercubic lattice 

z = Nd 1 E = ~vNd n(v) = Ns ^ N ) (53) 

here E is energy, and s(E) is entropy as function of the energy for corresponding ferromagnetic Ising model. At given 
0o w e can define corresponding v as minus energy per bond for ferromagnetic model at temperature -J-. 

vo0o) = (54) 
We can remember from thee definition of temperature 

< 55 » 

We obtain for the free energy 
0F 

-Jj2 = (1 - v 2 (0))Re\ncosh0 + s(v 2 (0)) + (v 2 (0) ~ v 1 (0))/2\ncosh20 1 (56) 

<- vlW , In 2 



+s(v 2 (0))/(2Nd)+0 1 dv nf(ir~J 
Jo 0(v o ) 



Integrating by parts in the last term we get 



0F 

^ = (1 - W2 (/3))i?elncosh/3 + s(v 2 {0))(v 2 {0) - « 1 (/3))/21ncosh2/3 1 (57) 



-\Ar-A-'» -.s(n(J.)))/(2.V,/)- C d0 o ^^+v 1 (0)0 1 y(0 1 ) 
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where y as a function of 0o is defined from the equation 



y = ff(^) (58) 
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function vq(/3) is denned from the (54). and j3\(0) is denned from the equation 



v 1 (/3) = -E0 1 )/(Nd) (59) 

where is defined from the (pl|). Here E0±) is the energy of ferromagnetic model on the same lattice. 

We propose the method for calculating approximate statistical mechanics of disordered models on finite lattices. It 
allows us also to identify Lee- Yang-Fisher zeroes. In the case of real temperatures our generalization of REM gave 
5% accuracy. If such an accuracy will be available here, then this method could be very efficient. 

We could succeed due to a simple observation (|4q) that on hierachy of generalized REM new (Lee- Yang-Fisher) 
phase could exist only on levels between paramagnetic and spin glass. On the other hand, there is no replica symmetry 
breaking. It is very interesting to consider the similar issues in strings where the ideas and methods of spin-glass 
theory have been applied recently. 
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